Spatial multipartite entanglement and localization of entanglement 
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We present a simple model together with its physical implementation which allows one to generate 
multipartite entanglement between several spatial modes of the electromagnetic field. It is based on 
parametric down-conversion with JV pairs of symmetrically-tilted plane waves serving as a pump. 
The characteristics of this spatial entanglement are investigated in the cases of zero as well as 
nonzero phase mismatch. Furthermore, the phenomenon of entanglement localization in just two 
spatial modes is studied in detail and results in an enhancement of the entanglement by a factor 

vN. 

PACS numbers: 03.67.Bg, 42.65.Lm 
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I. INTRODUCTION 

The existence of entanglement between more than two 
parties for discrete or continuous quantum variables is 
one of the most striking predictions of quantum me- 
chanics. The interest in multipartite entanglement is 
motivated not only by fundamental questions but also 
by its potential for applications in quantum communi- 
cation technologies. Possible applications of continuous- 
variable (CV) multipartite entanglement include quan- 
tum teleportation networks, quantum telecloning, con- 
trolled quantum dense coding, and quantum secret shar- 
ing (see Ref. [l| and references therein). 

Characterizing multipartite entanglement has raised 
much interest since the pioneering work of Coffman, 
Kundu and Wootters @. They have established for a 
three-qubit system and conjectured for ./V-qubit systems 
the so-called monogamy of quantum entanglement, con- 
straining the maximum entanglement between partitions 
of a multiparty system. More recently, Adesso, Serafini 
and Illuminati have introduced the continuous- 

variable tangle as a measure of multipartite entanglement 
for continuous-variable multimode Gaussian states. In 
particular, they have demonstrated that this tangle satis- 
fies the Coffman-Kundu- Wootters monogamy inequality. 
The conjecture of Ref. [2| has been proven by Osborne 
and Verstraete 0- The corresponding proof for Gaus- 
sian states has been obtained by Hiroshima, Adesso and 
Illuminati Q. 

A number of schemes to generate CV multipartite en- 
tanglement have been proposed theoretically and real- 
ized experimentally in recent years. There were first the 
passive optical scheme using squeezed states mixed with 
beam splitters @. Then came the active schemes in 
which multipartite entanglement is created as a result 
of parametric interaction of several optical waves such 
as cascaded/concurrent [Iol - fl2l |. interlinked [U, 13, or 
consecutive parametric interactions [l5j . These schemes 
generally neglect the spatial structure of the electromag- 
netic field. 

It is natural to investigate whether the spatial modes 
can also serve for the creation of CV multipartite en- 
tanglement. This question has been adressed recently in 



Ref. |16j | where a simple active scheme was proposed for 
the creation of tripartite entanglement between spatial 
modes of the electromagnetic field in the process of para- 
metric interaction. It consists in pumping a nonlinear 
parametric medium by a coherent combination of sev- 
eral tilted plane monochromatic waves which is called a 
spatially-structured pump. Since the pump photon can 
be extracted from any of these waves and the pair of 
down-converted photons emitted in different directions 
according to the phase-matching condition, that scheme 
allows for the creation of tripartite entanglement between 
spatial modes of the down-converted field. The aim of 
this paper is to generalize the results of Ref. [16[ to the 
generation of spatial multipartite entanglement. For that 
purpose, we present a scheme with an arbitrary number 
2N of tilted plane waves pumping a parametric medium. 
An interesting feature of this realistic proposal is the pos- 
sibility of localizing the created spatial multipartite en- 
tanglement in just two well-defined spatial modes formed 
as a linear combination of all the modes participating in 
the down-conversion process. The possibility of entan- 
glement localization was introduced by Serafini, Adesso 
and Illuminati [4j and a physical implementation was also 
given in terms of 2N-1 beam splitters and 2N single-mode 
squeezed inputs based on (l7j . 



The paper is organized as follows. Section II is devoted 
to spatial multipartite entanglement. We present the 
process of parametric down-conversion with a spatially 
structured pump consisting of N pairs of symmetrically- 
tilted plane waves in Subsec. II A. The evolution equa- 
tions for the field operators are given and solved in the 
rotating wave approximation for possibly non-zero con- 
stant phase mismatch. The genuine multipartite entan- 
glement is then studied In Subsec. II B. Section III is 
dedicated to the phenomenon of entanglement localiza- 
tion. It is presented explicitly in Subsec. Ill A and a 
quantitative characterization of the spatial entanglement 
generated by this process is given in Subsec. Ill B. The 
scaling with N as well as the dependence on the phase 
mismatch are explicitly provided. Conclusions are drawn 
in Section IV. 
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FIG. 1: (Color online) Scheme for the generation of spa- 
tial multipartite entanglement with 2N symmetrically tilted 
pump waves (illustrated for N = 2 and 4). The little circles 
represent the projections of the pump wave vectors in the xy 
plane of the crystal entering face. 



of the N pairs of opposed pumps on Fig. 1, for example 
the yz plane. 



— The evolution of the creation and annihilation opera- 

x tors, a and a\ associated to the propagation and diffrac- 
tion of the quantized electromagnetic field through the 
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II. SPATIAL MULTIPARTITE 
ENTANGLEMENT 

A. Parametric down-conversion with 
spatially-structured pump 

We consider a system of pumps which consists of N 
pairs of symmetrically tilted plane waves: 



E P (r) 



ft 



N 



(1) 



with r = (x,y), the vector of coordinates in the plane of 
the crystal entering face, and q = (k x , k y ), the projection 
of the three-dimensional wave vector k in that plane. Its 
Fourier transform in the xy plane of the (infinite) crystal 
(entering face) corresponds to Dirac delta's centered in 
q d et -q d 



N 



E p (q) 



£[*(q- 



qd) + S(q + q d ) 



(2) 



d=l 



This scheme is illustrated on Fig. 1 for N — 2 and 4. 
The little circles represent the projections of the pump 
wave vectors in the xy plane of the crystal entering face. 
The plane depicted in Fig. 2 is the one going through one 



FIG. 2: (Color online) The plane depicted is the one passing 
through one of the N pairs of opposed pumps on Fig. 1 (e. g. 
the yz plane). The tilted pumps have wave vectors k p (±q n ) 
with n ranging from 1 to N. The transverse (vertical) com- 
ponents are ±q n . 



nonlinear parametric medium is described by 18 1 

JU(z,q) = A J dq%(q-q')at(z,-q') e 4A ^'^. 

Here z is the longitudinal coordinate and A is the phase 
mismatch defined as 

A( q , - q ') = k z (a) + M-q') - Mq - q'). ( 4 ) 

where k z {q) and k z (—q) are the longitudinal component 
of the two incoming wave vectors and k pz is the longi- 
tudinal component of the pump wave vector. The con- 
servation of energy and momentum imply that we have 
to consider a set of 2 N + 1 coupled equations, obtained 
upon substitution of ^ into ([3]), for the waves exiting 
the crystal with wave vector transverse components q, 
<ld — q and — q^ — q with d = 1, • • • , N 
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dz 



N 

i(z,q) = aA^{a t (z,q d -q)e 4A(cl - cld - c ' )2 +a t (z,-q d -q)e 4A(q ^ qd - cl)2: } 



(5) 



d=l 
N 



dz 



a(z,±qj+q) = aA^ {^(z, q d T qj - q)e lA(q ' qdTqj " q) ^ + a f (z, -q d T qj - q)e iA(q '" qdTqj " q)2 } . 



d=l 



In the second equation, the phase mismatch is higher for 
the contributions with d ^ j which shall therefore be ne- 
glected. This corresponds to the usual rotating wave ap- 
proximation. The other phase mismatches are all taken 
equal to A which amounts to imposing some symetries 
on A(±qd ± q, ±q). One may introduce a renormalized 
phase mismatch 5 and a new relevant variable f = \f2a\z 
which combines the interaction strength aX and the lon- 
gitudinal coordinate z: 



A 



2 v / 2aA 
V2a\z. 



We shall use the following notation 



ao(r) 



a(z,q) 
a(z,±q„ 



(6) 



(7) 



71 = 1, 



,N, 



and shall consider only the zeroth spatial Fourier com- 
ponents of the held, q = 0. Physically, this corresponds 
to photodetection of the light field by a single large pho- 
todetector without spatial resolution. These modes are 
depicted in Fig. 2: the long (blue) arrows pertains to 
a n± (0), the short one to aj(0). In this setting, one can 
then rewrite ([5]) as 



d 
ar 



e 2iSf N 



d_ . 
df 



„2iSf 



(4+ + 4_) 



V2 



,N. 



(8) 



We may solve this set of equations and obtain for the 
fields at the output of the crystal, r = r\ z =e' 



ao(r) 
a n± (r) 



U(r)a o (0) 



V(r) 



a n± (0) 
V(r) 



/2N 
U(r) - 1 



JV 

E{4 + (o)+4jo)} 



d=l 

N 



'2N 



^2{a d+ (0) + a d _(0)} 



'2N 



4(o) 



d=l 

n = !,••• ,N. 



The functions U(r) and V(r), which satisfy |£/(r)| 2 
|y(r)| 2 = 1, are given by 



U{r) = e lSr ^cosh(\/iV7r) - i-^ sinh(VWjr)j 

V(r) = e i5r -sinh(\/7V7r), (10) 

7 



where 7 depends on the reduced phase mismatch S, 



N' 



(11) 



In the next subsection we shall characterize these solu- 
tions for the fields at the output of the crystal. 



B. Genuine multipartite entanglement 



The parametric down-conversion process preserves the 
Gaussian character of incoming modes. Hence, the out- 
going modes are Gaussian states which are thus com- 
pletely characterized by the covariance matrix associated 
to their quadrature components: 



Po(r) 
x n± (r) 
Pn ± (r) 



2Re a (r) 
21m ao(r) 
2Rea n± (r) 
21m a n± (r). 



7i= ,N 



(12) 



(9) From the solution (O one obtains 
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N 



t> j j I \" <a\ t TTt \- tn\ , X^Jr> r// n ^+ (°) + ^- (°) , T „, Jd+ (0) + Pd_ (0) 
Re (7(t-)xo(0) — Im U[r)po (0) + ^ RcV(r) — 1 ^== h Im K(r) 



d=l 

iV 



d=l 



M r ) 

Mr) = ImJ7(r)x (0) + ReU(r)p (0) +^jlmF( 
Pn± (r) 



V27V 

z d+ (0) (0) 
\/27V 



ReV(r\ 



V27V 



Re F(r) 
V27V 

Im V(r) 



Mo) 
Mo) 



JmV(r) 
ReV(r) 



Po(0) +x n± (0) 

M0)+Pn ± (0) 



N 



d=l 
N 



£N(Retf(r)-l 



^rf + (0) +xd_(0) 
27V 



/27V 



Im C/(r) 



E i Im£/(r) " d+(0) + ^- (0) + (ReU(r) 1) 



d=l 



Pd + (0)+Pd-(0) 
27V 

p d+ (o)+PdAo) 



27V 



27V 



(13) 



We can now determine the covariance matrix ele- 
ments ay = ^fA&A£,- + AfjAfi) /2\ of the out- 
put state p with (.) = Tr [p.] and A£j = ^ — 
£i ) where & is some component of the vector £ = 



(x ,po,ii + ,Pi + ,£i_,Pi_,£2 + ,P2 + ,--- ,x N _,p N _). Tak- 
ing into account that all inputs are vacuum states we 
obtain from (fT3")) 

([Mr)} 2 ) = ([po(r)} 2 )=a 
<[*n ± (r)] 2 ) = (^ ± (r)] 2 )=l 



a - 1 
27V 



(i n± (r)x d± (r)) = (i n± (r)x dT (r)) = 
<Pn±(r)pd ± (r)) = (p n± (r)p dT (r)) = 
(x (r)x n± (r)) = - (M r )pn ± (r)) = ■ 
(M.r)p n± ( r )) = (Po(r)x n± (r)) 



1 



2N 
a- 1 

27V 
b 

v/27V 
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(^o(^)Po(O) = (M( r )pd±(r)) = (x n± (r)M( r )) = 0, 

(14) 

for n, d = 1, • • • , TV, with 

a= |[/| 2 + |y| 2 

6 ee 2(Re U Re 7 - Im U Im V) 

c ee 2(Re U Im + Im URc V). (15) 

Ordering the lines and columns according to Xq, po, 

the covari- 



Pi+, xi_, pi_, x 2+ , P2+, 
ance matrix reads then 

/ A D D D 

D B C C 

D C B C 

D C C B 



D C C C 
D C C C 
\D C C C 



xn - , Pn- , 
D D D\ 

c c c 
c c c 
c c c 

B C C 
C B C 
C C B ) 



(16) 



where A, B, C and D are the following 2x2 matrices 



A = ah 



D 



B = 1 



1 



27V 



/27V 



h = 



C 

1 
1 



1 



27V 



(17) 

Its symplectic eigenvalues {z^} are \/ a 2 — b 2 — c 2 (with 
a fourfold degeneracy) and 1 (with a 27V — 3 degen- 



eracy). Noting that b 



= A\U\ 2 \V\ 2 implies that 



b 2 



1. The product of symplectic eigenval- 



ues is thus unity and, as expected for a unitary evolution, 
the output state remains pure. 

The covariance matrix (1161) is bisymmetric, i. e. 
it is invariant under the permutation of quadratures 
x n± (r),p n± (r) <H- x d± (r),p d± (£) or x d ^_ (r), p d ^ (£). As a 
consequence, it has the multipartite entanglement struc- 
ture of covariance matrices associated to bisymmetric 
(l+2N)-mode Gaussian states considered in Ref. [4j. We 
now show that the output state obtained here exhibits 
gen uine multipartite entanglement in the sense of Ref. 
(l7| . For that purpose, we have to verify that the follow- 
ing condition on covariance matrix elements is violated 



Q 




N 



E {*»+(**) +XnA r )} 
n=l 

N 

E {^«+( r ) +Pn-(r)} 




1 

27V' 
(18) 



From (fT4"|) - p5[) one deduces that 

Q = 4{a-b) (19) 
= 4 ([Re U(r) - Re V(r)} 2 + [Im U(r) + Im V(r)} 2 ) . 

When the phase matching condition is satisfied, 5 = 0, 
one deduces from (fTUl) that U(r) = cosh-\/7Vr and V(r) = 
sinh-\/7Vy. It follows that (|2ll|) reduces to 



Q = 4{cosh(V7Vr) - sinh(VTVr)} = 4e 



(20) 



5 



This quantity is smaller than 1/2N if the squeezing pa- 
rameter r satisfies 



= V2a\£ > 



31n(2AT) 



2VN 



(21) 



Under this condition on the pump amplitude a, the cou- 
pling parameter A and the crystal length £, the out- 
put state p produced by the above parametric down- 
conversion process with 2N symmetrically-tilted plane 
waves therefore exhibits genuine multipartite entangle- 
ment. This generalizes the result obtained in Ref. JT6J in 
the case of tripartite entanglement (N — 1). 

For a small nonzero phase mismatch, we can expand 
(121)1) around S — 0, which yields 

Q = Ae- 2 ^ r + ^ ([3 - 4iVr 2 ] e - 2v ^ 

+ 4[2ViVr - 1] + - l] 2 e 2VWr ^j + 0{S 4 ). 

(22) 

Owing to the last term, this quantity increases exponen- 
tially with y/Nr. As a consequence, Q remains smaller 
than 1/2N only for low values of the phase mismatch. 
This suggests that the genuine multipartite entangle- 
ment, at least when it is estimated with the criterion 
(|18[) is very sensitive to S. This situation is in contrast 
with the phenomenon that we study in the next section, 
namely the localization of entanglement which will be 
shown to be robust with respect to the phase mismatch. 



III. LOCALIZATION OF ENTANGLEMENT 
A. Beamsplitting the output state 

It has been shown Q that the entanglement of bisym- 
metric (to + n)-mode Gaussian states is unitarily local- 
izable, i. e., that, through local unitary operations, it 
may be fully concentrated in a single pair of modes. We 
shall study explicitly this phenomenon here. For that 
purpose we shall perform the following unitary transfor- 
mation based on discrete Fourier series, 



ao{f) 




-7zik 



: a n _ (f)} 
(23) 



Physically, this corresponds to beamsplitting the quan- 
tized fields at the output of the crystal. As a consequence, 
from (fT3"j) one deduces that the new fields at the output 
of the crystal are transformed to the following ones 



S'o(r) 
a[(r) 
a' k (r) 



U(r)a' (0) + V(r)a'S(0) 
U{r)a[(p) + V(r)atf{0) 
a' k (0) fe = 2,...,2AT, 



(24) 



Accordingly, the quadratures are now given by 



x' {r) = ReU(r)x' Q (0) ~lmU(r)p' (0) + ReV(r)x[{0) + lmV(r)p' 1 (0) 

p' (r) = ReU(r)p' o (0) +lmU(r)x' o (0) -ReV(r)p[(0) +lmV{r)x 1 (0) 

x[(r) = RcU(r)x[(0) - lmU(r)p[(0) + ReV(r)x' {0) + liaV(r)p' o (0) 

p[(r) = ReU(r)p[(0) +lmU(r)x' 1 (0) - Re V(r)#,(0) +lmV{r)x' o (0) (25) 

x' k (r) = 4(0) fc = 2,--. ,2N 

P' k {r) = p' k (Q). 



Hence, the covariance matrix elements associated to the 
density matrix p' are 



(Wo(r)} 2 ) = ([p'o(r)} 2 ) = <[4M] 2 ) = ([p[(r)} 2 ) = a 

(x> Q (r)x[(r))=-(p> (r)pi 1 (r))=b 

(x'o(r)Pi(r)) = (p'o(r)&i(r)) = c 

(Wk(r)] 2 ) = (Wk(r)] 2 )=l fc = 2,-..,2iV, (26) 



and the other ones are zero. Ordering the lines and 



p' 2 jy, the covariance matrix reads 

b c 



(a 

b 





a c 
c a 
-b 



2(2JV-1)) 



(27) 



columns according to x' , p' Q , x[, p[, x' 2 , p' 2 , 



l 2JV' 



where Ik is the k x k unity matrix. Its symplectic eigen- 
values {l/j} are the same as those of a as a' is obtained 
by congruence, 

a' = S T aS, (28) 

where the elements of the symplectic transformation S 
are given by (1231) . In the next subsection, we investigate 
the non-separability properties of the pertaining state p' . 



6 



B. Logarithmic negativity 

The covariance matrix (1271) is bisymmetric, i. e. the 
local exchange of any pairs of modes within its two diag- 
onal blocks leaves the matrix invariant Notice that 
here the unitary transformation (1231) is such that the co- 
variant matrix a 1 is block-diagonal, and invariant under 
the exchange of modes and 1 . 

The separability of the state p' can be determined as 
in Ref. [J] by (i) considering a partition of the system in 
two subsystems and (ii) investigating the positivity of the 
partially transposed matrix p' obtained upon transposing 
the variables of only one of the two subsystems. The 
positivity of partial transposition (PPT) is a necessary 
condition for the separability of any bipartite quantum 
state [H|,[2(|. It is also sufficient for the separability of 



(1 + n)— mode Gaussian states [21|,[22j. The covariance 
matrix a' of the partially transposed state p' with respect 
to one subsystem is obtained [2l| by changing the signs 
of the quadratures p'j belonging to that subsystem. 

Here, because of the block-diagonal structure of a' we 
can restrict the analysis to the first block and consider 



the transposition of mode 1, i.e. change the sign of p[: 

(29) 

Its symplectic eigenvalues {^} are 1 and 



( a be 

b -c a I /2 ^ 
\c b a 



v' ± = a ± V& 2 + c 2 = (\U\±\V\) 2 



(30) 



The necessary and sufficient PPT condition for the sep- 
arability of the state p' amounts to having v 1 - > lVj. We 
can therefore focus on the smallest eigenvalue v'_. The 
extent to which this criterion is violated is measured by 
£j^/(p'), the logarithmic negativity of p' defined as the 
logarithm of the trace norm of p': 

Em{p') =ln||p'||i =m&K(Q,-\nv'_), (31) 

where v'_ is obtained exlicitly from (|3"0T) and (ITU1) . 



1 - 



5 2 



N cosh* 



cosh ( V~N 7r ) — sinh ( V~N 



(32) 



r 



We first consider the case of zero phase mismatch, 5 
0, for which (|3"Tj) reduces to 



(cosh(VNr) - sinh(ViVr)) 



-2\/iVr 



(33) 



The logarithmic negativity is thus positive and, further- 
more, scales quadratically with the number of modes, 



(34) 



£ M (p') = 2VNr 



Note that the covariance matrix (|27fl is block-diagonal 
since (fT5]) entails that a = cosh(2-\//VV), b — sinh(2vMVr) 
and c = 0. It is readily recognized that its nontrivial 
part pertains to two entangled modes with a squeezing 
parameter 2\Nr and the remainder to 2N — 1 modes in 
the vacuum state. This is one instance of the localization 
of entanglement as introduced in Ref. \4 \ . The result (|3"4")) 
generalizes the one obtained in Ref. [16 1 for a single pair 
of tilted pump (N = 1) and shows that the effective 
squeezing is enhanced by a factor y/~N. 

When the phase mismatch takes on some finite val- 
ues, the logarithmic negativity ([3"Tj) features the smallest 
symplectic eigenvalue (|30l) . It is instructive to obtain a 
more explicit expression for a small phase mismatch by 



expanding around the case 5 = 
£m{ P ') = 2VNr- 



^tanWiVr 
N \ VNr 



The correction is of second order in S and negative. It 
is a bounded function of TV whose magnitude decreases 
for large N. The positivity of the logarithmic negativ- 
ity implies that the central mode is entangled with the 
uniform superposition of the tilted modes. This entangle- 
ment localization only slowly decreases when the phase 
mismath 6 increases. This confers some robustness to 
the entanglement localization process and is somehow in 
contrast with the genuine multipartite entanglement (|18[) 
which is more sensitive to the phase mismatch. Recalling 
that the squeezing parameter r = \f2a\l, equation ([33)1 
also quantifies explicitly the entanglement localization in 
terms of the pump amplitude a, the coupling parameter 
A, the crystal length I and provides the scaling with N. 



IV. CONCLUSIONS 

We have presented a simple active optical scheme for 
the generation of spatial multipartite entanglement. It 
consists in using 2N symmetrically tilted plane waves 



7 



as pump modes in the spatially structured parametric 
down-conversion process taking place in a nonlinear crys- 
tal. We have found the analytical solution of the cor- 
responding model in the rotating wave approximation 
for arbitrary N and possibly nonzero phase mismatch. 
We have studied quantitatively the entanglement of the 
2N + 1 coupled modes obtained at the output of the crys- 
tal. When the phase matching condition is satisfied, the 
system exhibits genuine multipartite entanglement. It 
has also been shown to subsist for nonzero, albeit small, 
values of the phase mismatch. 

In addition, our scheme provides a realistic proposal for 
the experimental realization of entanglement localization. 
By mixing different spatial modes with beam splitters, we 
can localize the entanglement distributed initially among 



all the 2N + 1 spatial modes in only two well-defined 
modes, formed by the linear combinations of the initial 
ones. Interestingly, this entanglement localization results 
in an enhancement of the entanglement by a factor \/N. 
Moreover, this process has been shown to be robust with 
respect to the phase mismatch. 
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